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Abstract. We present preliminary results from the analysis of an RXTE obser- 
vation of Cyg X-l in the hard state. We show that the observed X-ray spectrum 
can be explained with a model for an accretion disk corona (ADC), in which 
a hot sphere is situated inside of a cold accretion disk (similar to an advection 
dominated model). ADC Models with a slab-geometry do not successfully fit the 
data. In addition to the spectral results we present the observed temporal prop- 
erties of Cyg X-l, i.e. the coherence-function and the time-lags, and discuss the 
constraints the temporal properties imply for the accretion geometry in Cyg X-l. 


INTRODUCTION 

Cyg X-l is one of the most firmly established persistent galactic black hole 
candidates (BHCs). Its X-ray spectrum in the hard state can be described by 
a power-law with a photon-index T w (1.5- 1.7), modified by an exponential 
cutoff with a folding energy £ fo id « 150 keV [1,2] and reprocessing features 
below 10 keV. This spectral form has generally been interpreted as being due 
to an accretion disk corona (ADC). 

As it has been previously pointed out, the standard “slab-like” ADC has 
problems explaining the high-energy radiation of BHCs, since 1. the strength 
of the observed reflection component is too small to be consistent with a slab 
geometry [2], 2. the observed temporal characteristics of Cyg X-l cannot be 
explained with a pure slab-like ADC [3], and, 3., numerical models show that 
the spectral parameters of Cyg X-l are inconsistent with those predicted for 
slab-like ADCs [4,5]. 

In this paper, we use ADC spectra computed using a non-linear Monte Carlo 
code to model the observed spectrum of Cyg X-l. In addition, we present the 
global temporal characteristics of Cyg X-l. The results are based on 10 ksec of 
data from Cyg X-l obtained with the Rossi X-ray Timing Explorer (RXTE) on 
October 23, 1996. The data were reduced with ftools 4.0, using version 2.6.1 
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FIGURE 1. Spectrum of the best-fit sphere+disk model, and the ratio between the data 
and the best-fit model. The jump between the low-energy and the high-energy spectra 
is due to the uncertainty in the absolute effective areas of the instruments. The relative 
normalization between the HEXTE clusters and PCA is 0.70 in this fit. 

of the PCA response matrix, version 1.5 of the PCA background-estimator 
program, and the HEXTE response matrices dated March 20, 1997. Due to 
the uncertainy in the response matrix we model the spectrum between 6 and 
200 keV only and add a systematic error of 2% to the PCA data. 


SPECTRAL FITTING 

Spectral fitting of the RXTE data shows that the spectrum can be well 
described with a model of the form KE~ r exp(-E/E{), where T = 1.4 and 
Ef = 150 keV, plus an additional black-body with a temperature of (1.0 ± 
0.1) keV. Since we are ignoring the spectrum below 6keV the constraint on 
the temperature of the soft-excess is fairly weak. There is evidence only for 
a weak reflection component in the spectrum, with the covering factor being 
<0.1Q/2tT. 

We have computed grids of spectra for self-consistent ADC models using 
a non-linear Monte Carlo scheme based on the code of Stern et al. [6] and 
further described in [4,5]. We have written interpolation routines that enable 
us to fit the observational data with the data reduction software XSPEC [7]. 
The data and subroutines are available upon request. 




Modeling with spectra from a standard slab-corona, where the accretion 
disk is embedded between two hot ADCs, does not result in an acceptable 
fit. Due to the large covering fraction of the cold disk, a large amount of the 
coronal radiation is reprocessed in the cold disk, resulting in a large flux of 
thermal radiation which can effectively Compton-cool the corona. Therefore, 
the temperature of the slab ADCs cannot get high enough to produce a hard 
spectrum. Slab-models also have too strong of a soft-excess [4]. In contrast, 
the sphere+disk model, in which a spherical corona is situated at the center 
of a cold accretion disk, is able to explain the observed spectrum. In Fig. 1 
we show the best-fit obtained with the sphere+disk model, as well as the ratio 
between the data and the best-fit. The seed optical depth of the corona is 
r e = 2.1 ± 0.1 and the temperature of the corona is T = 65.7 ± 3.3 keV. The 
deviation between the model and the data is always less than 2%, which is an 
indication that the model is fitting the data quite well. Our formal x? e d' va l ue 
for this fit is 1.0. Note that the x 2 * va lue is affected by the statistical errors 
associated with Monte Carlo models, where computational time constraints 
prohibit the generation of models with a very high signal to noise ratio. 

TEMPORAL PROPERTIES 

Although there is good evidence from the spectral observations that the 
accretion geometry in Cyg X-l differs from a simple slab, the results do not 
rule out other geometries in which the covering fraction of the cold material is 
smaller than 30%. In order to further constrain the geometry we must consider 
the temporal properties of the source. For example, since in Comptonization 
the hard photons emerging from the source will have, on average, undergone 
more scatterings than low energy photons, there is a time-lag between hard 
and soft X-rays. The exact form of this time-lag can, in principle, be used to 
further constrain the accretion geometry of the system. Only a geometry in 
which both, the spectral and the temporal data, can be explained, should be 
considered a valid candidate geometry for Cyg X-l [8,9]. 

Much of the information on the geometry can be obtained by interpreting 
the phase information of the observed light curve [8,9]. In Fig. 2 we display 
the coherence between the energy band above 14.9keVand the energy band 
below 5 keV. The coherence-function for the frequency / is defined as [8] 

_ I<C(/>I 2 

7U> <|S,(/)| J x|S 2 (/)| 2 > 

where C(f) is the cross-spectrum between the intrinsic signals 5,(/) = 
Xi(f) - Ni(f) where X t (f) is the observed signal and N t (f) is the Pois- 
son noise-component (“observational rms noise”). The measured coherence 
of unity from 0.01 Hz to 10 Hz indicates a remarkable stability in the timing 
properties of the signal over the whole spectrum — an indication that either 




FIGURE 2. Coherence between the hard and soft energies. 


there is a single global source for the observed fluctuations, or that there is a 
global coherent response from the Compton corona, or both [3]. 

Similar information is also contained in Fig. 3, where the time-lag between 
the same energy bands is shown. In regions where the coherence function 
is unity, i.e., 0.01 Hz to 10 Hz, the time-lag can be well determined. Above 
~10Hz the coherence is not preserved and, Poisson noise begins to dominate; 
therefore, the time-lag varies erratically. 


Summary 

We have, for the first time, fit the RXTE spectrum of Cyg X-l using 
self-consistent ADC models. We have shown that the slab-geometry does 
not explain the spectrum of the object, while the sphere+disk model, having 
a smaller covering factor, can describe the observed spectrum. In addition, we 
have presented the measured coherence function and time-lags, which can be 
used to further constrain the accretion geometry. We emphasize that only the 
direct application of the models to the observational data can show whether 
a given geometry is applicable, or not. 

We acknowledge helpful discussions with M. Begelman, D. Gruber, I. Kreykenbohm, 
K. Pottschmidt, Ch. Reynolds, and R. Staubert. This work has been financed by NSF 
grants AST91-20599, AST95-29175, INT95-13899, NASA Grant NAG5-2025, NAG5-3225, 
NAG5-3310, DARA grant 50 OR 92054, and by a travel grant to J.W. from the DAAD. 



Time Lag for (14.09-100 keV) vs. (0-3.86 keV) 



FIGURE 3. Time lag measured between the hard and the soft energy band. Stars indicate 
a negative time lag, i.e. the hard band leads the soft band. The minimum observed time 
lags of < 1ms indicates a maximum coronal size of < O(50GM/c 2 ) for M = 10 M©. The 
maximum observed time lag constrains the dynamics of the corona. 
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We present the results of the analysis of the broad-band spectrum of Cygnus X-l from 
3.0 to 200 keV, using data from a 10 ksec observation by the Rossi X-ray Timing Explorer. 
Although the spectrum can be well described phenomenologically by an exponentially 
cut-off power law (photon index T = 1.45+q!o 2 > e-folding energy Ef = 162*|keV, plus a 
deviation from a power law that formally can be modeled as a thermal blackbody with 
temperature kT BB = 1.21J;® keV), the inclusion of a reflection component does not improve 
the fit. As a physical description of this system, we apply the accretion disc corona models 
of Dove, Wilms & Begelman [4] — where the temperature of the corona is determined 
self- consistently. A spherical corona with a total optical depth r = 1.6 ± 0.1 and an 
average temperature kT c = 87 ± 5 keV, surrounded by an exterior cold disc, does provide 
a good description of the data (Xred — 1-55). These models deviate from the data by 
up to 7% in the 5-10 keV range, and we discuss possible reasons for these discrepancies. 
However, considering how successfully the spherical corona reproduces the 10-200 keV 
data, such “photon- starved” coronal geometries seem very promising for explaining the 
accretion processes of Cygnus X-l. 


1. Introduction 

Cygnus X-l is one of the most firmly 
established persistent galactic black 
hole candidates (BHCs). The combi- 
nation of a Comptonization continuum 
and reprocessing features has been in- 
terpreted as being due to an accre- 
tion disc corona (ADC). The geomet- 
ric configuration of the corona and 
the cold disc, however, is still un- 
clear. In most prior work, the ge- 
ometry has been assumed to be a 
cold accretion disc embedded between 
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two hot coronae, with a slab ge- 
ometry [10,9] .Recently, however, evi- 
dence has been presented showing that 
ADC models with a slab geometry 
suffer from several problems, making 
them less likely to be the appropriate 
models for explaining the high-energy 
radiation of BHCs. Gierlitiski et 
al. [8] analyzed simultaneous Ginga- 
OSSE data of Cyg X-l and found 
that the strength of the reflection com- 
ponent is too small to be consistent 
with a slab geometry. Using a non- 
simultaneous broad-band spectrum of 
Cyg X-l, we previously argued that 
our self-consistent ADC models with 
a slab geometry (in which the coronal 
temperature is not a free parameter 
but is determined by balancing Comp- 
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ton cooling with viscous dissipation) 
could not have a high enough tem- 
perature to explain both the observed 
power-law and exponential cutoff [4,5]. 

In this paper, we analyze the broad- 
band spectrum of Cyg X-l from 3.0 
to 200 keV, using data from the Rossi 
X-ray Timing Explorer (RXTE). The 
simultaneity of this broad-band ob- 
servation allows us to place stronger 
constraints on the coronal parameters 
than previously possible. In §2, we 
give the results from spectral modeling 
of the Cyg X-l spectrum using both 
“standard” simplified models as well as 
our self-consistent ADC models. In §3 
the results of these fits are discussed. 

2. Spectral Analysis 

RXTE observed Cyg X-l for a total 
of 22.5 ksec on 1996 October 23 and 
24. For the analysis presented here, we 
use data from the proportional counter 
array (PCA) and the high energy X- 
ray timing experiment (HEXTE). 

2.1. Standard Models 

Spectral fitting was first performed 
using the following standard models: 
a power-law, a power-law with an ex- 
ponential cutoff, a power-law with an 
exponential cutoff plus a cold reflec- 
tion component, and thermal Comp- 
tonization models. In addition, we 
added a Gaussian line (with energy 
and width fixed to 6.4 keV and 0.1 
keV, respectively) to several of these 
models. We fixed the low-energy ab- 
sorption to an equivalent cold Hydro- 
gen column of N H = 6 x 10 21 cm -2 , the 
value suggested by the soft X-ray spec- 
trum and interstellar reddening mea- 
surements of HDE 226868 [1,14]. For 
energies E 10 keV, a power-law with 
an exponential cutoff provides a good 
description of the data, as also shown 
in Figure 1. The value for the e-folding 
energy, Ef , is well constrained due to 
the quality of the HEXTE data. The 
analysis of the residuals of this fit in- 
dicates the presence of a soft excess at 
energies & 8keV, as shown in Figure 1. 
Adding a black-body component with 
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Figure 1. Comparison of the RXTE 
spectrum of Cyg X-l to the best-fit 
power-law with an exponential cutoff 
model. This model was fit to only the 
10-200 keV data. The photon power- 
law index T = 1.46 and the e-folding 
energy Ef = 167 keV. A soft-excess be- 
low lOkeV is clearly evident. 


a temperature kT BB « IkeV and a 
Gaussian component, with an equiva- 
lent width EW= 46.4 eV, to the model 
improves the quality of the fit, result- 
ing in a good description of the data 
over the whole energy range (Xr e( j = 
169/166), as shown in Figure 2. We 
note that this black-body component 
should not be interpreted literally as 
being due to disc emission: it is simply 
added phenomenologically to measure 
the magnitude of the data’s deviation 
from a power law. We could have in- 
stead used a disc blackbody, Gaussian, 
or bremsstrahlung model with roughly 
equal success. 

The addition of a reflection com- 
ponent to this model does not im- 
prove the quality of the fit, as the 
best-fit value for the fraction of in- 
cident radiation that is Compton re- 
flected by cold matter is / ^ 0.02. 
Even though we only considered PCA 
data between 3.0 keV and 30keV, the 
power-law index of the continuum is 
well constrained by the HEXTE data, 
and there is no evidence that the spec- 
tral slope of the data for the two in- 
struments disagree within the energy 
range of the overlap. The PCA residu- 
als between 20 and 30 keV seen in Fig- 
ure 1-4 are consistent with those seen 
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spectrum of Cyg X-l to the best-fit 
power-law with an exponential cutoff 
model. The photon power-law index 
T = 1.45 and the e-folding energy 
is Ef = 163 keV. A blackbody com- 
ponent, with a temperature /cXbb = 
1.1 keV and a relative flux of & 5%, 
and a Gaussian component, with an 
energy fixed at 6.4 keV and an equiv- 
alent width EW= 46.4 eV, were added 
to fit the soft-excess. For this fit, 
x? ed = 169/166. 


for power-law fits to the Crab. 

It is this 30-100 keV energy range 
that is crucial for constraining a 
power-law index since these energies 
are uncontaminated by reprocessing 
or reflection features. We therefore 
postulate that the observed 10-30 keV 
“hardening” feature typically found by 
others appears to be the beginning of 
a hard power-law that continues to 
~ 100 keV, and that the ~ 1-10 keV 
portion of the spectrum has a softer 
power-law due to contamination from 
the thermal excess. This contamina- 
tion could be due to Comptonization 
of the thermal radiation emitted by the 
cold disk with an effective blackbody 
temperature kT hB ~ 150 keV, since, for 
energies & 10fcr BB , the Comptonized 
spectrum is not expected to be a power 
law. 

We believe that the measured spec- 
trum parameters, which represent a 
harder power law and a weaker re- 
flection component than previously 
found by others, is mostly due to 
our fitting the entire 3-200 keV spec- 


tral band. Previous observations us- 
ing different instruments have had rel- 
atively poor spectral coverage between 
~ 30-100 keV. results are heavily in- 
fluenced by any remaining uncertain- 
ties in the PCA response matrix. To 
support this claim, if we restrict the 
data to 3-30 keV we do find best-fit 
parameter values similar to those pre- 
viously found for Cyg X-l. Specifi- 
cally, for this restricted data range, we 
find that a reflected power law with 
T = 1.81 ± 0.01 and covering fraction 
/ = 0.35 ±0.02 yields *? ed = 157/90 = 
1.75. Adding back the HEXTE data in 
the 30-100 keV range significantly re- 
duces the quality of the fit (x 2 e d = 8.7, 
77 dof). Therefore, being able to fit a 
reflected power-law with a large disk 
covering fraction is seen to be partly 
a consequence of ignoring the data in 
the ~ 30-100 keV range. We refer the 
reader to our forthcoming paper [6] for 
a more detailed discussion of this issue. 

2.2. Accretion Disc Corona mod- 
els 

For both a slab ADC model and 
the sphere+disc ADC model, we have 
computed grids of spectra using a non- 
linear Monte Carlo scheme based on 
the code of Stern et al. [13] and mod- 
ified by Dove et al. [4] and Dove 
et al. [5]. Slab-corona models do 
not result in good fits, as the pre- 
dicted spectra are always much softer 
than the observed spectrum (x 2 ed ~ 
80). This result is consistent with 
our previous findings based on non- 
simultaneous data [5]. 

The sphere+disk ADC model does 
provide a good description to the data 
(cf. table ?? and figure 3). The formal 
X 2 -value of our best-fit model (x 2 e d = 
1.55) is larger than the values found for 
some of the phenomenological models 
discussed in §3.1. Considering that our 
model has only three parameters (seed 
optical depth, coronal compactness, 
overall normalization, plus the “hid- 
den” parameter Xbb> which is fixed 
at 150eV during the fitting proce- 
dure), and considering that our model 
is physically self-consistent , the level of 
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Figure 3. Top: Best-fit sphere+disc fit 
to the PCA and HEXTE data. Bot- 
tom: Ratio between data and model 
for the best-fit sphere+disc model. 
Deviations less than 2% are within 
the range of uncertainty of the PCA 
response matrix. Deviations in the 
20-30 keV range are similar to those 
seen with best-fit models of Crab neb- 
ula observations (Jahoda 1997, private 
communication) . 


agreement with the data is quite good. 
The major disagreement between the 
model and the data occurs in the « 5- 
10 keV range, where the residuals are 
as large as 7%. 

3. Conclusions 

We have applied a variety of spec- 
tral models to an RXTE observation 
of Cyg X-1. We find that the observed 
spectrum between 3.0 and 200 keV is 
well described by a power-law with 
a photon index T = 145^o!o2 m °d- 
ified by an exponential tall with an 
e-folding energy Ef « 160 keV, and 
a soft excess that was modeled as a 
black-body (although other broad dis- 
tribution models are capable of ex- 
plaining the soft-excess) having a tem- 
perature kT BB « 1.2 keV. The mea- 
sured value of r is lower than those 
found in previous broad-band analy- 
ses, which find a T ~ 1.65 [8,3]. In 
addition, the observed strength of the 
Compton reflection feature is weak, as 


the best-fit covering fraction, deter- 
mined by fitting the data with a power- 
law reflection model [11], is found to 
be « 0.2 when no soft excess compo- 
nent is included, and is found to be 
& 0.02 when one is included. The lat- 
ter value is considerably lower than the 
/ « 0.3 found by Gierlinski et al. [8]. 
The sphere+disc ADC models provide 
a good explanation of the data, and 
furthermore these models are physi- 
cally self-consistent. For a more de- 
tailed discussion of these results, we 
refer the reader to Dove et al. [6]. 
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